Can inhomogeneities solve the horizon problem ? 
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Abstract 

We show how inhomogeneous cosmological models can naturally explain the large angle correlation 
we observe in the CMB (cosmological microwave background) radiation without invoking any infla- 
tionary stage, but simply considering the effects of inhomogeneities on the propagation of photons 
from the last scattering surface. 
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I. INTRODUCTION 



High redshi 
surement [3 



uminosity distance measurements 19|, |25|, |26|, [27|, [28|, [29] and WMAP mea- 
3l| of cosmic microwave background (CMB) interpreted in the framework of 
standard FLRW cosmological models have strongly disfavored a matter dominated universe, 
and strongly supported a dominant dark energy component, corresponding to a positive cos- 
mological acceleration, which we will denote as a FLRW (not to be confused with the scale factor 

n 

a). As an alternative to dark energy, it has been proposed [6] that we may be at the center 
of a inhomogeneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) solution 
of Einstein's fields equations, where spatial averaging over one expanding and one contracting 
region is producing a positive averaged acceleration ap, even if it was shown that quantities 
constructed in this way can be unobservable because they violate the causal structure of the 
underlying space. Another more general approach to map luminosity distance as a function 
of redshift Dl(z) to LTB models has been recently proposed [21], where an inversion method 
was applied to successfully reproduce the observed Dl(z), also showing that the freedom in the 
choice the LTB geometry poses some degeneracy problem, which could be solved by imposing 
other cosmological constraints such as CMB radiation, or RSSE (redshift spherical shell energy) 

3- 

The main point is that the luminosity distance is in general sensitive to the geometry of the 
space through which photons are propagating along light geodesies, and therefore arranging 
appropriately the geometry of a given cosmological model it is possible to reproduce a given 
Dl(z). For FLRW models this correspond to set constraints on and Q m and for LTB models 
it allows to determine the functions E(r), M(r),tf,(r). 

In a similar way the inhomogeneities determine a different causal structure from the ho- 
mogeneous case, in particular a spatial dependency of the Hubble radius H(r)~ l , which can 
allow causal contact of points at the surface of last scattering, without invoking a period of 
exponential expansion. 

A previous attempt to solve the horizon problem in LTB spaces was investigated in 
,but this approach involved a special choice of the bang function £&(r) , while in our case we 
consider a homogeneous big bang scenario, i.e. = 0. Since the time scale of observation 

of CMB is very short on cosmological scale, of the order of few years, its observed isotropy is 
compatible with a spherically symmetric radially inhomogeneous model, but it is not a sufficient 
evidence of homogeneity. 
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II. LEMAITRE-TOLMAN-BONDI (LTB) SOLUTION 

Lemaitre-Tolman-Bondi solution can be written as 



as: 



(R, r f dr 2 



ds 2 



-dt 2 + 



+ R 2 dtt 2 



(1) 



l + 2E{r) 

where R is a function of the time coordinate t and the radial coordinate r, E(r) is an arbitrary 
function of r, and R, r denotes the partial derivative of R with respect to r. 
Einstein's equations give: 

2 2E(r) 2M(r) 

' + -^r 1 , (2) 
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with M(r) being an arbitrary function of r and the dot denoting the partial derivative with 
respect to t. The solution of Eq. ([2]) can be written parametrically by using a variable rj = 
J dt/R, as follows 

M(r) 



R(r],r) 



-2E(r) 

M(r) 
-2E(r) 



1 - cos ( ^-2E(r)r] 

sin (j-2E(r)rj) 
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(4) 
(5) 



where R has been introduced to make clear the distinction between the two functions R(t, r) 
and R(rj, r) which are trivially related by 



R(t,r) = R(ri(t,r),r) (6) 

and tfi(r) is another arbitrary function of r, called bang function, which corresponds to the 
fact that big-bang/ crunches happen at different times in this space. This inhomogeneity of the 
location of the singularities is the origin of the possible causal separation between the central 
observer and the spatially averaged region for models with positive ay- 

We can introduce the following variables 

6M(r) 



. . R(t,r) . 2E(r) 
a(t,r) = -^, k(r) = -±1, p {r) 



(7) 



so that Eq. ([T]) and the Einstein equations (|2J) and Q can be written in a form which is more 
similar to FRW models 
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ds 2 = -dt 2 + a 2 



1 + 



a, r r 



a 



dr^ 



1 — k(r)r 2 
k(r) p (r) 



+ r 2 dQl 



+ 



3a 3 



p{t,r) 



(por 3 ), r 
6a 2 r 2 (ar) 



The solution in Eqs. (j3j) and (jSJ) can now be written as 

Po{r 



a(fj,r) 
t(fj,r) 

where fj = rjr = J dt/a . 



6k(r) 

Pop) 
6k(r) 



1 — cos ( \J k(r) 7] 
1 

rj - 



k(r) 



sin I J k(r) fj 



+ U{r) 



(9) 
(10) 

(11) 
(12) 



III. PROPAGATION OF PHOTONS AFTER LAST SCATTERING IN A INHOMO- 
GENEOUS SPACE 

When we study cosmology in a not homogeneous space it is important to adopt the appro- 
priate distance and time redshift relations and use them consistently. In particular the time 
of last scattering tf s RLW and tj^ rB , are different , but in the following we will omit the upper 
script for the cosmological model, assuming we are referring to a LTB model. 

We will assume that the Universe was homogeneous at the time of last scattering and that 
inhomogeneities affected only the propagation of photons from the last scattering surface until 
us today. In order to compute U s and r; s is is necessary to integrate backward in time the null 
geodesic equation as a function of redshift. The geodesic equations expressed in terms of the 



redshift can be written as 



dr 



l + 2E(r(z)) 



dz (1 + z)d t d r R(t(z),r(z)) 
dt -\d r R(t(z),r(z))\ 



(13) 
(14) 



dz (1 + z)d t d r R(t(z),r(z)y 
where t(z) and r(z) physically represent the coordinates along the null geodesic of the photon 
coming to us (located at r = 0, z = 0). After integration we get: 



ris = r(z ls ) 

Us — t(zi s ) 



(15) 
(16) 
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which in general will be different from those obtained in a FRLW space. 
Points which are today at an angular distance a c were at the time of last scattering ti s at a 
distance: 



d ts = a c R(t ls ,r ls ) (17) 

If we focus on a constant r spatial hyper surface, we can see that the notion of event horizon 
is the same as in a homogeneous space, but it depends on r: 



v ; a R(t,r) V 



As long as 



d l8 < H-\n 8 ) (19) 
pw^j- 1 (20) 

the correlation of CMB on angular scale a c will compatible with the causal structure of 
inhomogeneous space. 

A stronger constraints would come by imposing that the Hubble horizon is less than the 
sound horizon, but we will limit to impose the upper limit given by ([20]) . 



IV. FLAT UNIVERSE 

We will now study the case of a flat LTB space, i.e. when E(r) = . Einstein equations 
can be solved analytically: 

R(t,r) = (^2M(r)[t-t b (r)]) 2/3 (21) 

Imposing the homogeneous bang condition tb(r) = 0, eq. (!20|) gives the following constraint on 
M: 

M(r ls ) < || (22) 

The energy density in this case is 

which corresponds to a FRW homogeneous matter dominated universe. This upper limit 
is in agreement with the general case discussed in the next section, which requires M(r) to 



be a decreasing function of r as a general feature of LTB models which can solve the Horizon 
problem. 



V. CURVATURE DOMINATED UNIVERSE 



The other limiting case which can be treated exactly is the curvature dominated Universe 
in which M(r) = . Also in this case it is possible to find an exact solution for R: 

R(t,r) = y/2E(r)[t-t b (r)) (24) 

The condition ( 120|) in this case is: 

E{n s ) < ^ (25) 

This upper limit is in agreement with the general case discussed in the next section, which 
requires E(r) to be a decreasing function of r as a general feature of LTB models which can 
solve the Horizon problem. 



VI. GENERAL CASE 

Using the implicit solution (jlj) we can find a general expression for the condition (1201) . in 
terms of rj. 

Defining a(r) : 



/ dR(r],r)\ 1 dt(r],r) 1 

\ dV ) d V J-2E(r ls ) 



tan 



rn s J-E(r ls )/2 



(26) 



we obtain : 



a{ri s ) > a c 



(27) 



For any given LTB model there will be a different rju corresponding to the time of last 
scattering, which can be found by solving numerically the equation ((SJ) 



Us = t(r)i s ,ri s ) (28) 

The quantity a(r) is the present angular projection of the Horizon scale H[r, t(r)] -1 , where 
t(r) is the time along a null geodesic which arrive at r = at present. In other words it is 
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the maximum present angular scale which could be in causal contact at co-moving coordinate 
r and time t(r). The r dependence is the crucial feature of LTB spaces, and it is what can be 
used to solve the Horizon problem with an appropriate choice of E(r) and M{r). 

The general inequality (|27|) is apparently independent of M(r), but it in fact depends im- 
plicitly through r]i s , since t(rj,r) contains M(r). For E(r) = and M(r) = equation fl28|) is 
indeterminate, so the general condition (1271) cannot be applied to these limit cases, which justify 
our previous separate treatment. 

As it can be seen in figure (El) the solution of eq. (128j) is inversely proportional to M(r), and 
a(r) is proportional to r]i s and inversely proportional to E(r). This implies that decreasing 
E(r) and M(r) are the right type of functional dependence on the comving coordinate to solve 
Horizon problem. 
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FIG. 3: In order to show the dependence of rji s on LTB geometry the solution of eq. ([28j) is plotted for 
t is = 10 , 0.3 < M(r) < 10 and 0.3 < E(r) < 5. 

In order to provide further insight on the dependence of the a(r) and its dependence on r)i s 
and E(r) we plotted it for different cases. 

The crucial feature of a(r) is that it is a decreasing function of E{r) , and its upper limit 
is r]i s /2 that is reached as E(r) goes to zero, which implies that in order to solve the Horizon 
problem E(r) has to be a decreasing function of r. 

As an example, fig.® shows how the choice E(r) = 1/r can solve the Horizon problem 
assuming r]i s = 2tt, if r/ s is sufficiently large. From this example it is clear that the important 
parameters to determine if a given model explains the observed large angle correlation of CMB 
radiation are r ts , ti s , r] [s , which depend implicitly on the the LTB geometry through the geodesic 
equations (TH)) . This pose in general a non trivial numerical problem, but our qualitative an- 
alytical analysis provides the main feature that the functions E(r), M(r) should satisfy, and 
which could be used as general directions in a more systematic numerical investigation. 

VII. CONCLUSION 

We have shown that inhomogeneous cosmological models can in principle not only be alter- 
natives to dark energy but also provide a potential alternative solution to the horizon problem, 
and we have set up the constraints that solutions have to satisfy. 

Both dark energy and inflation correspond in the framework of FRLW cosmology to a stage 
of positive acceleration, and it consequently possible that if homogeneities can mimic dark 
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energy they could also provide the causal structure necessary to solve the Horizon problem, 
which could then be considered an artifact originating from neglecting inhomogeneities effects 
on the propagation of photons from the last scattering surface. 

It must be observed though, that there may not be any model satisfying the different cosmo- 
logical constraints and the conditions we have derived, and in general the two sets of conditions 
could be incompatible. Despite that, it is important to observe that inhomogeneities could play 
an important to explain any physical situation which treated in a FRLW framework involves a 
stage of positive acceleration. 

The extensive search for a model which successfully solves the horizon problem and is com- 
patible with other cosmological observables will be the subject of a future investigation. 
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